library(bpca)
## Loading required package: scatterplot3d
## Loading required package: rgl
## Warning: package 'rgl' was built under R version 3.2.3
library(rgl)
library(ggfortify)
## Warning: package 'ggfortify' was built under R version 3.2.3
## Loading required package: proto
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 3.2.3
library(ggplot2)
library(RColorBrewer)
library(knitr)
library(rglwidget)
knit_hooks$set(webgl = hook_webgl)

setwd("/Users/bobminnich/Documents/Columbia/Courses/Data_Mining/Examples/DigitReader")
data = as.data.frame(read.csv("train.csv", header = TRUE, sep = ","))
labels2 = data[,1]
labels_frame <- as.data.frame(data[,1])
labels_frame <- setNames(labels_frame, c("l"))
data_r = data[,-1]
pca = prcomp(data_r)
plot.new()
screeplot(pca, main = "PCA Plot of NIST data", type = "lines")

dev.off()
## null device 
##           1
colorpal = c("#E41A1C", "#0066ff", "#4DAF4A", "#984EA3", "#FF7F00", "#FFFF33", "#A65628","#ff37cb","#66ff33", "#00ffff")

#Find colors associtated with labels and apply the color palet
for(i in 1:10){
  labels_frame$color[data$label == i-1] = colorpal[i]
}
#Used {r testgl, webgl=TRUE, } for R Chunk
#Plotting
plot3d(pca$x[,1:3],col = labels_frame$color, size = 1)
#Set initial View
rgl.viewpoint(theta = -35, phi = 30, fov = 45, zoom = 1)
#Add Legend (Only viewable in R)
legend3d("topright", legend = paste('Num', c(0:9)), pch = 16, col = colorpal, cex=1, inset=c(0.15))

#Legend 
pie(rep(1,10),labels = c(0:9), col=colorpal, main = "Hand Written Digits - Legend")
Legend : Scroll to Zoom, Click and drag to rotate
Num 0 Num 1 Num 2 Num 3 Num 4
Num 5 Num 6 Num 7 Num 8 Num 9

You must enable Javascript to view this page properly.